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Abstract 

We investigate the scahng of Tsalhs entropy in disordered quantum spin-S chains. We show that an extensive 
scahng occurs for specific values of the entropic index. Those values depend only on the magnitude S of the 
spins, being directly related with the effective central charge associated with the model. 
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1. Introduction 

Correlations among parts of a quantum system 
are behind remarkable phenomena, such as a quan- 
tum phase transition (QPT) [H In particular, 
the relationship between correlations and QPTs is 
revealed by the behavior of entanglement at criti- 
cality as measured, e.g., by the von Neumann en- 
tropy (see, for instance, Ref. Q). Given a quantum 
system in a pure state \ip) and a bipartition of the 
system into two subsystems A and B, the von Neu- 
mann entropy between A and B reads 



S = -Tr {pA In pa) = -Tr {ps Inps) , 



(1) 



where pA — Tr^p and pb — Tr^p denote the re- 
duced density matrices of A and B, respectively, 
with p — HA and B are probabilistic in- 

dependent (such that p = PA ® Pb), the von Neu- 
mann entropy is additive, i.e.. Sab = Sa + Sb- 
As a consequence, S is extensive for uncorrelated 
subsystems, namely, S{L) cx L, where L denotes 
the size of a block of the system. On the other 
hand, S becomes nonextensive in presence of corre- 
lations. Indeed, for critical systems in one dimen- 
sion, which are known to be highly entangled, con- 
formal invariance implies a diverging logarithmic 
scaling given by S{L) cx (c/3) InL (or, more specif- 
ically, S{L) — (c/3) InL -I- constant), where c is the 
central charge associated with the Virasoro algebra 
of the underlying conformal field theory 0, [1, 
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For noncritical (gapful) systems in one dimension, 
entanglement saturates at a constant value k, i.e., 
5(L)^>fcasL— >oo. More generally, for higher di- 
mensions, noncritical systems are expected to obey 
the area law, which implies that the von Neumann 
entropy of a region scales as the surface area of 
the region instead of the volume of the region it- 
self. In other words, the area law establishes that 
S{L) cx L^^^ (L CX)), where D is the dimension 
of the system. 

Remarkably, it has recently been shown in 
Refs. [7, 8] that a quantum system may exhibit 
specific probability correlations among its parts 
such that an extensive entropy can be achieved even 
for highly correlated subsystems. This has been ob- 
tained by generalizing the von Neumann entropy 
into the nonadditive Tsallis q-entropy 



1 



(Trp* - 1) , 



(2) 



with q G R. One can show that the von Neumann 
entropy is a particular case of Eq. ^ by taking 
(7 = 1. Tsallis entropy has been successfully ap- 
plied to handle a variety of physical systems, in 
particular those exhibiting long-range interactions. 
Recent experimental results for its predictions can 
be found, e.g., in Refs. ll|, [l^. In Tsallis statis- 



tics, the additivity of the von Neumann entropy for 
independent subsystems is replaced by the pseudo- 
additivity relation of the Sq entropy 



Sq[pA "S) Pb] = SqlpA] + Sq[pB] 
+ (1 - q) Sq[pA]Sq[pB]■ 



m 
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The investigation of Sq in conformal invariant quan- 
tum systems has revealed that the extensivity of the 
entropy can be achieved for a particular choice qext 
of the entropic index q in Eq. ([2]). In particular, 
Qext is directly associated with the central charge c. 
More specifically, the extensivity of Sq occurs for 0] 



Qext 



V9 



(4) 



The aim of this work is to consider the scaling of 
the nonadditive entropy Sq and, consequently, its 
extensivity in quantum critical spin chains under 
the effect of disorder into the exchange couplings 
among the spins. Indeed, disorder appears as an 
essential feature in a number of condensed matter 
systems, motivating a great deal of theoretical and 
experimental research (see, e.g., Refs. 13) 3]) ■ 
particular, it is well known that, in the case of a 
spin-S random exchange Heisenberg antiferromag- 
netic chain (REHAC), disorder can drive the sys- 
tem to the so-called random singlet phase (RSP), 
which is a gapless phase described b y sp in singlets 
distributed over arbitrary distances [l5|. In recent 
years, it has been observed that the entanglement 
entropy in critical random spin chains displays a 
logarithmic scaling that closely resembles the be- 
havior of pure (non-disordered) systems. Indeed, 
for a block of spins of length L, we have that the 
von Neumann entropy reads S{L) cx (ce///3)lnL, 
where Ceff is an effective central charge that gov- 
erns the scale of the entropy 16]. Moreover, it has 
been shown that in the case of the RSP, Ce// is 
determined solely in terms of the magnitude S of 
the spin in the chain [13, [H, 111 (see Ref. 2^ for 
a review of entanglement in random systems and 
Ref. 2l| for other connected results). Here, we will 
show that the extensivity of Sq can also be obtained 
for random critical spin chains, with q^xt governed 
by Cef f ■ Hence, Qext will be given as a unique func- 
tion of the spin S. Moreover, as we will see, around 
the extensivity point qext, Sq{L) cx L'^^ with the 
exponent 7 of the power law given by a quadratic 
function of q. 



Figure 1: A schematic picture of the RSP. Spin singlets are 
composed randomly at arbitrary distances. 



In order to evaluate Sq in the RSP, we begin by 
considering a number n of singlets connecting a con- 
tiguous block composed by L spins with the rest of 
the chain. In this situation, the pseudo-additivity 
of Sq implies that Tsallis entropy is given by the 
Proposition below. 

Proposition 1. For a bipartite system composed 
of a number n of spin-S singlets connecting two 
blocks, with n G N, Tsallis entropy Sq^^ for each 
block is given by 
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(25+ l)"(i-«) - 1) 



(5) 



Proof. The proof can be obtained by finite in- 
duction. Indeed, the single-site reduced density 
operator pA for a spin-S singlet can be repre- 
sented by a D-dimcnsional diagonal matrix given 
by PA = diag{D-^,D-^,---,D-^), with D = 
25-1-1. Therefore, from Eq. (0), we obtain that 
= (1 - q)-^iD^-i - 1). For two singlets, the 
pseudo-additivity of Sq given by Eq. ([3]) implies 
that 5f ^ = (1 - q)-\D^'^^-i^ ~ 1). By taking the 
general expression for the entropy for n singlets as 
S^q^ = (l-g)-i(i:)"(i-9)-l), we obtain for (n-hl) 
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singlets that 5^"+'^ = (1 - g)-i(Z?("+i)(i-9) 
Hence, Eq. ([5]) holds for any 71 e N. ■ 

Tsallis entropy for the RSP can then be ob- 
tained by numerically averaging iSq"'' over a sam- 
ple of random couplings along the chain. These 
random configurations are generated by following 
a gapless probability distribution, which drives the 
system to the RSP, with the entropy of each config- 
uration computed by counting the spin singlets via 
a renormalization group approach described in the 
next section. 



2. Nonadditive entropy for a set of random 
singlets 

We begin by considering the typical arrange of a 
quantum spin-S chain in the RSP, which is provided 
by a set of spin singlets distributed over arbitrary 
distances, as sketched by Fig. [TJ 



3. Renormalization group method for ran- 
dom spin systems 

The RSP can be conveniently handled via a per- 
turbative real-space renormalization group method 
introduced by Ma, Dasgupta and Hu (MDH) 0, 
[23l |. which was successfully applied to the spin- 1/2 
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REHAC. This approach was proven to be asymp- 
toticahy exact, which ahowed for a fully character- 
ization of the properties of the RSP Q^. Consid- 
ering a set of random Heisenberg antiferromagnetic 
interactions Ji between neighbouring spins 5*^ and 
Si+i, the original MDH method consists in finding 
the strongest interaction D, between a pair of spins 
{S2 and S3 in Fig. 2a) and treating the couplings of 
this pair with its neighbors ( Ji and J2 in Fig. 2a) 
as a perturbation. Diagonalization of the strongest 
bond leads at zeroth order in perturbation theory 
to a singlet state between the spins coupled by n. 
Then, the singlet is decimated away and an effec- 
tive interaction J' is perturbatively evaluated. By 
iteratively applying this procedure, the low-energy 
behavior of the ground state will be given by a col- 
lection of singlet pairs and the structure of the RSP 
will naturally appear. 
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the strongest coupled pair obtaining an effective in- 
teraction between the neighbors to this pair which 
is given by Eq. ^ (see Fig. 2a). This new effective 
interaction is always smaller than those eliminated. 
Now suppose Ji > J2 and Ji > 3n/[2S{S + 1)]. In 
this case, we consider the trio of spins-S* coupled by 
the two strongest interactions of the trio, Ji and fl 
and solve it exactly (see Fig. 2b) . This trio of spins 
is then substituted by one effective spin interacting 
with its neighbors through new renormalized inter- 
actions obtained by degenerate perturbation theory 
for the ground state of the trio. This method has 
been used to successfully investigate the quantum 
phase diagram of the spin-1 ^2^ and spin-3/2 EgI ] 
Heisenberg spin chains and will be here applied to 
the computation of the nonadditive entropy. 

4. Scaling of nonadditive entropy in the RSP 

Let us apply now the generalized MDH approach 
to analyze the behavior of Sq in the RSP of antifer- 
romagnetic spin chains. We start with a REHAC, 
whose Hamiltonian is given by 



N 



(7) 



Figure 2: Modified MDH renormalization procedure for spin- 
S chains. 

Unfortunately, when generalized to higher spins, 
this method, at least in its simplest version, re- 
vealed to be ineffective. The reason is that, af- 
ter the elimination procedure of the strongest bond 
O, the effective interaction J' may be greater than 
ri. Then, the problem becomes essentially non- 
perturbative for arbitrary distributions of exchange 
interactions. For instance, considering the REHAC 
with arbitrary spin-S', the renormalized coupling is 
given by the recursive relation p3 | 



J' 



1) 



J1J2 



(6) 



Notice that, for S* > 1, the renormalization fac- 
tor is (2/3)S'(S' -I- 1) > 1, resulting in the break- 
down of perturbation theory. In order to solve this 
problem, a generalization of the MDH method was 
proposed in Refs. [2^, [2^ (for other proposals, see 



also 27 



28|). This generalized MDH method con- 
sists in either of the following procedures shown in 
Fig. [5] Taking the case of the Heisenberg chain as 
an example, if the largest neighboring interaction to 
fi, say Ji, is Ji < 3ri/[2S'(S'+l)], then we eliminate 



where {Ji} are random exchange couplings obey- 
ing a probability distribution P{J) and {"i^i} are 
spin-S operators, with periodic boundary condi- 
tions adopted. The numerical investigation of Sq is 
performed as follows. We begin by considering spin 
chains with 200, 000 sites, whose couplings {Ji} are 
randomly generated by using a gapless power law 
distribution P{J) oc J""'*, for which trio renormal- 
izations are negligible. Results are then obtained by 
averaging over a sample of M = 40, 000 configura- 
tions for {Ji}0- For each random configuration j, 
we decimate the spins out via the generalized MDH 
technique and compute the number Uj of singlets 
that cross a block of len gth L 
counted by following Ref. 



The number Uj is 
17[ . In turn, if a singlet is 



decimated and the spins composing the singlet are 
in different blocks, this singlet adds one to the total 



number 



On the other hand, in the case of a trio 



elimination, nothing is added to nj, since one effeo 
five spin is returned to the chain (see also Ref. [l^ 



^It is worth observing that a considerably larger amount 
of configurations is demanded for the evaluation of Tsallis 
entropy in comparison with the von Neumann entropy (see, 
e.g., Ref. 0). 
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Figure 3: (Color online) Nonadditive entropy Sq{L) as a 
function of the block size L for a spin- 1/2 REHAC for q = 
-1.30 until q = -1.50. 
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Figure 4: (Color online) Nonadditive entropy Sq{L) as a 
function of the block size L for a spin-1 REHAC for q = 
-0.40 until q = -0.60. 



for a similar approach). The decimation procedure 
is iterated until the elimination of all spins in the 
chain. After all configurations computed, Tsallis 
entropy is then obtained by averaging over all of 
them, i.e., 



(8) 



with iSq"^^ given by Eq. The results for spin- 
1/2, spin-1, and spin-3/2 REHACs are shown in 
Figs. [3][5l Note that an extensive Sq can be 
found for specific negative values Qext in the ex- 
amples above. These values are clearly distinct for 
different spin magnitudes S. In turn, a remarkable 
fact is that q^xt depends only on S for a system in 
the RSP and not on the specific model considered. 



Figure 5: (Color online) Nonadditive entropy Sq{L) as a 
function of the block size L for a spin-3/2 REHAC for q = 
-0.1 until q = -0.3. 



In fact, let us consider the spin-1 random exchange 
biquadratic antiferromagnetic chain, whose Hamil- 
tonian is given by 



H 



Biq 



N 



(9) 



Application of the generalized MDH procedure here 
results only in the formation of singlets, with the 
renormalized exchange coupling reading (23 | 



J'- 



2 J1J2 
9 Vl ' 



(10) 



where Ji and J2 are the nearest neighbors of the 
strongest bond Vi. As in the Heisenberg case, we 
consider spin chains with 200, 000 sites, whose cou- 
plings {Ji\ are randomly generated by using a gap- 
less power law distribution P{J) oc J~°'* and aver- 
aged over a sample of M = 40,000 configurations 
for { Ji}. This model produces the same result for 
qext as the spin-1 REHAC in the RSP. In order to 
illustrate this equivalence and obtain a numerical 
value for qext^ let us consider the relationship be- 
tween Sq and L close to the extensivity index qext- 
As indicated in Figs. [3][5l we obtain that Sq oc L'^ 
around the extensive point 7 — 1. Indeed, the aim 
of Figs. [3J5] is precisely show this dependence as q 
is varied across extensivity. Moreover, as shown in 
Fig- m for the REHACs (with different spin magni- 
tudes S) as well as for the spin-1 biquadratic chain, 
the exponent 7 shows a quadratic dependence on 5, 
namely. 



^ — uq + vq + w, 



(11) 
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Figure 6: (Color online) Exponent 7 of the power law 
Sq{L) oc as a function of q. We fit 7 as a quadratic 
function of g, with the extensivity index q^xt varying only 
with the magnitude S of spin in the chain. In particular, note 
the collapse of the spin-1 REHAC and biquadratic curves, 
indicating that q,,xt is uniquely determined by S. 

with u,v,w e R. Note in Fig. [5] the coUapse of 
the spin-1 curves for different models, which indi- 
cates that Qext is a unique function of the spin of 
the chain. Therefore, we can associate qext with the 
effective central charge that governs the entangle- 
ment scaling in the RSP, which is given by [l3| 

Ce//=ln(25 + l). (12) 

Moreover, q^xt can be directly obtained by impos- 
ing 7 = 1 in Eq. (jlip . The values for q^xt as well as 
their relationship with the effective central charge 
are summarized in Table [T] Concerning the esti- 
mation of the error bar Aq^xt, it can be directly 
obtained from a standard error propagation proce- 
dure. Indeed, from A7 = \dj/dq\ Aq, we obtain 
Aqext = A"f / \ 2uqext + v\. Although the analytical 



Model 


qext 


Ceff 


Spin-1/2 REHAC 


-1.40 ±0.03 


In 2 


Spin-1 REHAC 


-0.49 ±0.04 


In 3 


Spin-1 biquadratic 


-0.55 ±0.06 


In 3 


Spin-3/2 REHAC 


-0.25 ±0.02 


In 4 



Table 1: Entropic indices q^xt that yield the extensivity of Sq 
for different magnitudes S and their corresponding effective 
central charges. 



relation between qext and S (and therefore Ce//) 
remains elusive at the present stage, it is reason- 
able to suppose that qext monotonically increases 
to unity (Boltzmann-Gibbs-von Neumann entropy) 
as S increases to infinity, where the classical limit 
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Figure 7: (Color online) Extensivity index qext as a function 
of the effective central charge c^f j for REEIACs in the RSP. 

is obtained. This behavior is indeed displayed in 
Fig. El in terms of 1/ce// for REHACs in the RSP. 
In particular, note that a linear behavior can be 
inferred between qext and 1/ce//, which reads 

1.67 

qext = 1 • (13) 

Therefore, as given by Eq. @ for the pure case, we 
can also determine an expression for qext in terms 
of Ceff for disordered systems, which reinforces the 
universal properties of the entropic index q in crit- 
ical systems. 

5. Conclusions 

In summary, we have investigated the scaling of 
Tsallis entropy Sq in spin-S random critical quan- 
tum spin chains. By focusing on the RSP, we have 
shown that, for specific values qext of the entropic 
index, Sq becomes an extensive quantity, which rec- 
onciles the quantum scaling with the Clausius-like 
prescription for classical thermodynamics. It is im- 
portant to emphasize that the extensivity of the 
nonadditive entropy does not imply absence of cor- 
relations between the parts of the quantum sys- 
tem. In this context, there is no contradiction with 
the behavior of block entanglement as measured by 
the von Neumann entropy. Remarkably, qext is di- 
rectly associated with the effective central charge 
Ceff that governs bipartite entanglement in ran- 
dom spin chains, which means that qext is solely 
determined by the magnitude S of the spin in the 
chain. Moreover, we have inferred a linear algebraic 
relationship between qext and 1/ce//. An analyti- 
cal investigation of this relationship as well as the 
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behavior of q^xt in other critical phases of random 
spin chains (e.g., Griffiths phase) are challenges un- 
der research. Moreover, the scaling of Sq far from 
the extensivity point is also relevant in connection 
with the scaling of Trp* (see, e.g., Refs. 29l 30j). 



We intend to address such topics in a future work. 
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